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In this note, we extend an evolutionary stochastic portfoho optimization framework to 
include probabilistic constraints. Both the stochastic programming-based modeling envi- 
ronment as well as the evolutionary optimization environment are ideally suited for an 
integration of various types of probabilistic constraints. We show an approach on how to 
integrate these constraints. Numerical results using recent financial data substantiate the 
applicability of the presented approach. 



1 Introduction 



o 
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Stochastic programming is a powerful method to solve optimization problems under uncertainty, 
see Ruszczyiiski and Shapiro 2003 for theoretical properties and Wallace and Ziemba 2005 



for an overview of possible applications. One specific application area is portfolio optimization, 
which was pioneered by H.M. Markowitz Markowitz 1952 . The advantage of using the stochastic 



programming approach is that the optimization can be done without using a covariance matrix 
of the assets, which is on one hard to estimate and on the other hand does not capture the 
uncertainty in sufficient detail, especially if there are many assets. 

Instead of using the asset means and the covariance matrix, a stochastic programming ap- 
proach uses a set of scenarios, each weighted by a certain probability. In the specific portfolio 
optimization context one scenario is a set of one possible asset return per asset under consider- 
ation - see below for more details or e.g. see Chapter 16 of Cornuejols and Tiitiincii 2007 . 



This note is organized as follows: Section [2] summarizes the evolutionary approach, which 
was used to solve the stochastic portfolio optimization problems. Section |3] adds probabilistic 
constraints to the standard problem and shows an approach on how to integrate these type of 
constraints Section [4] summarizes numerical results using recent financial data, while Section [5] 
concludes the note. 



2 Evolutionary Stochastic Portfolio Optimization 



We follow the approach taken by |Hochreiter 2007 and Hochreiter] [2008| , which builds a stochas- 
tic programming-based environment on top of the general evolutionary portfolio optimization 



findings reported by Streichert et al. 2003 , Streichert et al. 2004a , and Streichert et al. 2004b . 
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2.1 Stochastic portfolio optimization 



Let us define the stochastic portfoho optimization problem as follows. We consider a set of assets 
(or asset categories) A with cardinality a. Furthermore, we base our decision on a scenario set 
S, which contains a finite number s of scenarios each containing one uncertain return per asset. 
Each scenario is equipped with a non-negative probability Ps, such that X^sgS^** ~ "^^^ 
scenario set can be composed of historical data or might be the output of a scenario simulator. 

For every portfolio x we can easily calculate the profit and loss distribution by multiplying 
the scenario matrix with the portfolio weighted by the respective probability. Let us denote the 
profit function of a certain portfolio x by Tr{x) and the loss function by £{x) — —tt{x). 

Every portfolio optimization procedure is a multi-objective optimization. In the traditional 
case it is a trade-off between return and risk of the profit and loss function. We do not want 
to employ an multi-objective approach such that we use the classical Markowitz approach and 
use the Variance of the loss distribution for the risk dimension, which we want to minimize, and 
the expectation of the profit function for the reward dimension, on which we want to set a lower 
limit. Hence, the main optimization problem is shown in Eq. ([I]) below. 

minimize Variance (i'a;) 
subject to E{tTx) > M 

Furthermore we consider the classical constraints, i.e. budget normalization, as well as setting 
an upper and lower limit on each asset position, as shown in Eq. ([2|. These are naturally fulfilled 
by the evolutionary approach shown below, especially since we restrict short-selling in our case. 

subject to Eae.A^a = 1 ,'2^ 
I <Xa <u yae A ^ ' 



2.2 Evolutionary stochastic portfolio optimization 

The evolutionary algorithm chosen is based on the commonly agreed standard as surveyed by 



Blum and Roli [2003' and is based on the literature cited at the beginning of this Section. 

We use the following genetic encoding of a portfolio. Each gene consists of of two parts: 
One that determines the amount of budget to be distributed to each selected asset and one part 
which determines in which assets to invest. The first part gi consists of a predefined number b 
of real values between and 1 and the second part 172 is encoded as a bit-string of the size of the 
amount of assets. 

The following evolutionary operators have been implemented and used: 

• Elitist selection. A certain number oi of the best chromosomes will be used within the 
next population. 

• Crossover. A number 02 of crossovers will be added, both 1-point crossovers {gi and 32) 
and intermediate crossovers (only gi). 

• Mutation. 03 mutated chromosomes will be added. 

• Random addition. Furthermore 04 randomly sampled chromosomes are added, which are 
also used for creating the initial population. 

The specific number of operators o = (01,02,03,04) has to be determined for each concrete 
number of assets a as well as the parameter b. 
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3 Probabilistic Constraints 



The main advantage of the stochastic programming approach is that the complete distribution 
is naturally available and can be used for optimization purposes by integrating these directly 
into the constraint handling mechanism. In the most simplest case, we want to restrict that the 
probability that the loss exceeds a certain value S is lower than a specified probability e. Given 
our profit function n^, the constraint we want to add to our optimization model is given in Eq. 

subject to F{nx < S) < e. (3) 

We will not treat probabilistic constraints as a hard constraint, but use the size of the violation 
for adding a penalty to the objective function. Let the fitness value which we aim to minimizing 
be /. We calculate the penalty p by 

P^fx (P(7r, < ^) _ e) X 7, 

where 7 is a factor to control the penalization level. The fitness value used for evolutionary 
optimization purposes is thus given by /' = / + p. Such a constraint can be implemented and 
handled conveniently. 



4 Numerical Results 

The program code was implemented using MatLab 2008b without using any further toolboxes. 

We use data from the 30 stocks of the Dow Jones Industrial Average at the beginning of 
2010, i.e. the ticker symbols AA, AXP, BA, BAG, GAT, GSGO, GVX, DD, DIS, GE, HD, HPQ, 
IBM, INTG, JNJ, JPM, KFT, KG, MGD, MMM, MRK, MSFT, FEE, PG, T, TRV, UTX, VZ, 
WMT, XOM. Daily data from each trading day in 2009 has been used. Weekly returns have 
been calculated. 

We are using b = 100 buckets which are distributed to the respective asset picks, such that 
each chromosome has a length of 6 + a = 130. The initial population consists of 1000 random 



chromosomes. The operator structure defined in Section 2.2 is o = (100, 420, 210, 100), such that 
each follow-up population has a size of 830. This number is due to the different combinations 
between the crossovers and mutations of gi and 32- 

Eirst, we optimize without using the probabilistic constraints, i.e. the main optimization 
problem given Eq. [l] using /i — 0.001. Then we add the probabilistic constraint shown in 
Eq. |3]with 5 = and e — 0.1. These values have been chosen after analysing the resulting loss 
distribution. 

Eig. [1] shows the optimal portfolio without applying the probabilistic constraint Pi (left) and 
the optimal one using the probabilistic constraint P2 (right). The allocation changed consider- 
ably. The diversification has not changed, i.e. three assets (KET, VZ, XOM) are dropped from 
the portfolio, and three others are picked (HPQ, MSFT, MMM). The resulting loss distributions 
are shown in Eig. [2] where the impact of the probabilistic constraint is immediately visible. 
Furthermore, the statistical properties of the portfolios are shown in Table [T] In this table, the 
naive 1/A^ portfolio P3 has been added for comparison purposes. While P3 has the highest re- 
turn, it is also the most risky one in terms of both risk parameters - standard deviation and the 
probability to fall below the specified threshold. Another interesting fact is that the probabilistic 
constrained portfolio yields a higher expected return than the standard optimal portfolio. This 
is partly due to the fact that the lower level /i has been set to a level below the expected return 
of the standard portfolio but is definitely another indicator that the plain classical Markowitz 
approach should not be used for contemporary portfolio optimization purposes. 
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Figure 1: Pi without P constraint (left) and P2 with P constraint (right). 




Figure 2: Loss distribution of Pi (left) and P2 (right). 



Pi (no P) P2 P3{l/N) 

Mean 0.0024 0.0051 0.0062 

Std.Dev. 0.02 0.0225 0.0398 

Prob. 0.1774 0.1089 0.2702 



Table 1: Statistical properties of various portfolios. 
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5 Conclusion 



In this note, an extension of an Evolutionary Stochastic PortfoHo Optimization to inchide prob- 
abihstic constraints has been presented. It can be shown that the integration of such constraints 
is straightforward as the underlying probability space is the main object considered for the op- 
timization evaluation. Numerical results visualized the impact of using such constraints in the 
area of financial engineering. Future extensions of this work include the integration of risk mea- 
sures into the probabilistic constraint, e.g. constraining the maximum draw-down of the optimal 
portfolio. 
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